Electromagnetic excitation rates for nuclear isomers in a hot 

dense plasma 



T.C Luu, J.L. Friar, and A.C Hayes 
Theoretical Division, Los Alamos National Laboratory 
MS-227 Los Alamos, NM 87545, USA 



Correspondence: Thomas Luu 

Los Alamos National Laboratory 
T-6 Division, MS B227 
Los Alamos, NM 87545 
Phone: (505) 667-3612 
Fax: (505) 664-0007 
tluu@lanl.gov 



18 total pages 
3 tables 
1 figure 



1 



LA-UR-04-5803 



Electromagnetic excitation rates for nuclear isomers in a hot 

dense plasma 

T.C.LuuQ J.L. Friar, and A.C. Hayes 
Theoretical Division, Los Alamos National Laboratory, 
MS-227, Los Alamos, New Mexico 87545, USA 
(Dated: February 7, 2008) 

Abstract 

In high neutron flux environments where isomers can be strongly populated by nucleonic re- 
actions, isotope abundances from reaction network chains can be affected by the population of 
nuclear isomers. At high temperatures and densities there is the additional possibility of populat- 
ing these isomers electromagnetically. Here we examine the rates for electromagnetic excitation of 
the isotopes of several isomers of interest both in astrophysics and applied physics (e.g. ^'^^U, ^^'^Ir, 
and s'i'jSSY), \Yg consider six possible electromagnetic processes, namely, photo-absorption, inverse 
internal conversion, inelastic electron scattering. Coulomb excitation, (7, 7') and (e, e'7) reactions. 
We find that for plasma temperatures kT ~ 1 — 10 keV the electromagnetic reactions rates are 
negligible. Thus, we conclude that reaction network calculations do not need to include for the 
possibility of electromagnetically exciting nuclear isomers. This is true in both stellar and terres- 
trial thermonuclear explosions, as well as in plasma conditions expected at the National Ignition 
Facility. 
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I. INTRODUCTION 



In a hot dense plasma short-lived excited states of the nucleus with excitation energies 
on the order of the temperature of the plasma reach thermal equilibrium with the nuclear 
ground state and their relative population is determined by a Boltzmann distribution. Nu- 
cleosynthesis calculations usually try to take these states into account, though little is known 
about their nucleonic reaction properties (see, for example, Ref[l]). Recently there has been 
increased interest in the reaction properties of a certain set of nuclear excited states: nuclear 
spin isomers. These isomers are characterized by their relative long lifetimes due to their 
large angular momentum compared to their respective ground states. Particular interest 
has been concentrated on the 77 eV isomer of ^^^U, where studies suggest that this isomer 
is populated strongly in the (n, n') reaction as well as in the ^^^U(?t,, 7) reaction, and that 
the fission and neutron capture cross sections for this isomer are significantly different than 
that of the ^^^u ground stateQ, 3- Thus, in high neutron flux environments, where two- or 
multi-neutron reactions are highly probable, reactions on this isomer can affect the reaction 
network chain involving the Uranium isotopes. The question we wish to address here is 
whether nuclear isomers can be populated by electromagnetic interactions as well. In a hot 
dense plasma environment there are several electromagnetic mechanisms to be considered 
at leading and sub-leading order in the fine structure constant a, namely: photo-absorption, 
inverse internal conversion, inelastic electron scattering. Coulomb excitation, (7,7'); and 
(e,e'7)Q]. Figure n shows the Feynman diagrams for these processes. 

The electromagnetic excitation rates for all of these processes depend, of course, on the 
temperature kT of the plasma, the excitation energy of the isomer, and the lifetime 
Ti/2 of the isomer. To simplify calculations, we assume that nuclei are completely stripped 
of their electrons in the plasma. Where possible, we give estimates to the errors of this 
assumption. We present expressions for these electromagnetic processes and examine the 
rates for certain nuclear isomers. Of particular interest for applied physics are the low-lying 
isomers of ^^^U, ^^^Ir, '^'^'^^Y. Their properties, as well as those of other isomers considered 
in this paper, are listed in Table B The half-lives of these isomers range from 0.3 msec to 
years and their excitation energies from ~ 77 eV - 2.5 MeV. 
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II. ELECTROMAGNETIC EXCITATION RATES 



The total reaction width can be expressed as the sum of its constituent widths: 



i.e.s. 



(1) 



Traditionally, this width is expressed as 




total 



r^(l + ai,c. + ai,e.s. + ...)• 



(2) 



Here i.e. refers to internal conversion, i.e.s. refers to inelastic electron scattering, and so 
forth. As can be seen from Eq. calculation of a particular channel width Fj involves 
calculating both and a^, where = 1 for the case of photo- absorption. The reaction rate 
per nucleon for a particular channel i is given by 



where is the flux of incoming particles responsible for that particular reaction and ai{E) 
is the cross section. Due to the relatively long lifetimes of these isomers {i.e. T total ^ -E'm), 
the cross section is well approximated by 



where Em is the excitation energy of the isomer and k is the wavenumber for the incoming 
particle. 

A. Photo-absorption 

We assume that in the system of interest the plasma has reached thermal equilibrium so 
that the photon flux is described by a Plankian distribution. 




(3) 





(4) 
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eN{E)dE, 
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where N{E) is the photon number density in the energy window between E and E + dE. 
The photo-absorption rate per unit atom is then 



R^ = c N{E)cT^{E)dE (6) 
Jo 

2 r 

= tsL (r) 

Here cr^{E) is the photo-absorption cross section. Equation [7| follows from substituting Eq. E] 
into Eq. El In the long-wavelength approximation (which is valid for the plasma conditions 
considered here) the width can be expressed in terms of reduced multipolarity transition 
probabilities B^ttX), 

In general, calculating B{7cX) is very difficult because it depends on the initial- and final-state 
nuclear wavef unctions, which are not typically known. Naive estimates can be made if one 
invokes the extreme single-particle approximation, giving the well-documented Weisskopf 
expressions jsl, 

2 



B{MX : J, ^ J,) = H (^)'i?^^-^ /.^fm^ 



(9) 



where i? ~ 1.2 A^^^ and /xat is the Bohr magneton. Near closed-shell nuclei, Weisskopf 
estimates give reasonable agreement with experimental values. However, for nuclei that 
exhibit strong collective degrees of freedom, the Weisskopf estimates are too small, sometimes 
by a few orders of magnitude 0]. 

Since the half-life of the isomer of ^35^ jg about 26 minutes, temperatures of ~ 150 
keV would be required for the reaction rate per unit atom to reach a value of one-per- 
second. Thus, photo-absorption does not provide a significant mechanism for populating 
the 77 eV isomer of ^^^U. As can be seen from Table IIIIl the photo-absorption rates for all 
the isomers considered here are negligible. Even with more accurate Weisskopf values, the 
photo-absorption rates of these isomers at these temperatures would still be very small. 
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B. Inverse Internal Conversion 



In inverse internal conversion an electron from the plasma is captured into a bound atomic 
orbital, thereby exciting the nucleus. The capture energy of the electron must match the 
excitation energy of the nuclear isomer. The reaction rate is determined by the overlap of 
the continuum electron wavefunction with the bound atomic electron wavef unction, weighted 
by the nuclear transition matrix element. Since we assume that nuclei are stripped of all 
their electrons, all atomic shells {i.e. K, Lj, Lu, Mj, etc) can participate in the internal 
conversion process, as long as the kinematics of the process are satisfied. However, since 
the conversion happens near the nucleus, s-waves states typically dominate the transition. 
In most cases the electron can capture into the lsi/2 -ft'-shell. For isomers with very low 
excitation energies the electron can be captured only into a higher orbit. For example, the 
77 eV isomer of 235U can capture electrons only into very high (n > 39) outer atomic orbits, 
whereas both ^^^Am and ^^'^Ir can capture electrons only into orbits with n >2. 

To obtain naive estimates of a^.c, we assume relativistic hydrogenic wavefunctions for the 
bound electron and treat the nucleus as a point particle. The internal conversion coefficients 
for the innermost atomic orbits are thus given byjj] 

2 

|2 



2 2 



2 



(10) 



a... (MA; '^o = -1) = X:(2^- + !)(«: - 1)M J ^ | |i?-i,.(MA)p, 



2 



where j = — 1/2 and is the photon wave number. For MX transitions, k, can take on 
values of —A or A + 1, while for EX transitions, k, = X or —(A + 1). The functions Rko,k are 
given by 

R.UEX) =iK, - k){Rs + R4) + A(i?4 - R3) + XiRi + R2) ^^^^ 
i?«„,,(MA) =R, + Re, 



6 



where 



oo 



Jo 

POO 

R2= dr v^{per)V^^{Xer)h^^\u^r) 
Jo 

POD 

i?3 = / dr v^{per)U^^{Xer)h^^}_^{uj^r) 

■^}o. ^ (12) 

Jo 

POO 

i?5 = / dr Vf,{per)Uf,^{Xer)h'-^\uJ^r) 
Jo 

POO 

i?6 = / dr u^{per)V,,^{Xer)h''^\u^r). 
Jo 

Here U and V correspond to the large and small components of the relativistic radial wave- 
function for the bound electron in a Coulombic field, while u and v are the corresponding 
analogues for the Coulomb-distorted continuum wavefunction, and hx is the spherical Han- 
kel function of order A. Closed analytic forms can be derived for these integrals and are 
given, for example, in Refs. P, 0] . 

Internal conversion coefficients for the innermost s-wave reactions are given in Table ITTl It 
is interesting to note that though internal conversion is suppressed by an order a compared 
to photo-absorption, there can be instances where there are enhancements to the reaction 
rate. This is clearly evident in ^^^Ir, ^^^Am, and especially ^^^U. This behavior, though not 
apparent from Eq. EH comes from the observation that the coefficients scale roughly asjlflj] 

A+5/2 



a,; 



(^) for£;A 



1 



A+3/2 

for MX 



Since the excitation energies of these isomers are relatively small, this enhancement is not 
surprising. 

In the systems of interest {e.g. ignited inertial confinement fusion (ICF) capsules, nu- 
clear explosions, and stellar environments), the temperatures attainable by the plasma are 
typically not high enough for nuclei to be completely stripped of their electrons, so that 
not all atomic shells are accessible for inverse internal conversion. Screening effects aris- 
ing from the partial stripping of these nuclei affect the isomer excitation reaction rates by 



as much as 50% or morellTtl. Self-consistent Dirac-Hartree-Fock calculations 



ll| have been 



performed to obtain internal conversion coefficients on some of the isomers considered here. 



These calculations include the full effects of screening since the nuclei are not considered 
stripped of their electrons. Finite-size effects of the nucleus are considered as well. Hence 
these calculations use more realistic electronic wavefunctions than the simplistic hydrogenic 
wavefunctions considered here. We list these results in Table ITTl for comparison. In general, 
calculations can differ by an order of magnitude (or more), as is evident for the case of ^^^Y. 
Yet the overall trend in relative sizes is consistent between the two methods. 
The density of electrons is given by 



rif. 



V27c^h^ J ^^e(^-^)/'=^ + 1 ^^^^ 

Since there is no simple closed expression for the integral above, the chemical potential n 
and density must be numerically solved self-consistently. The flux $e of electrons is simply 



47^2^3 J e(E-M)/fcT + l 

Hence the differential flux is 



f ,^ E ^ ^ 

dE- 



d(^,) = dE (15) 



Using Eqs. OH gives the rate for isomeric production via internal conversion as 



(16) 



4 ^{Em-tl)/kT + 1 ■ 

To obtain estimates, we take an average electron density of rif. ~ IQQNa cm~'^, where 
A'^ is Avogadro's number. This corresponds to a chemical potential of /i = —.00275 keV 
(-.0622 keV) at fcT = 1 keV (10 keV). Rates are listed in Table |Th1 Even with the more 



accurate internal conversion coefficients of Ref. 



ll|, rates would still be neghgible. 



C. Inelastic electron scattering 

Formally, there is little difference between the theory of inelastic electron scattering and 
internal conversion. Physically, the main difference is that the final state electron still 
belongs to the continuum, as opposed to the bound-state orbital characteristic of inverse 
internal conversion. If the point-nucleus approximation is made, calculations of 'inelastic 
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scattering coefficients' are identical to the case with internal conversion coefficients-one just 
replaces the bound state wavefunctions with Coulomb distorted continuum wavef unctions, 
i.e. 



Ri = 


/ dr u^{pr)u^^{por)h^p 
Jo 


{u^r) 


R2 = 


POO 

/ dr Vf,{pr)Vf,XPor)h^x'< 
Jo 


[u^r) 


R3 = 


POO 

/ dr v^{pr)u^XPor)h^x- 
Jo 




i?4 = 


POO 

1 dr u^{pr)v^^{por)hfl 
Jo 






POO 

/ dr v^{pr)ut,XPor)h^x 
Jo 


{u^r) 


Re = 


POO 

/ dr u^,{pr)v^,XPor)h^x 
Jo 


{u^r) 



(17) 



The above integrals can be computed numerically and the coefficients calculated using 
Eq. E3 Note that since these calculations deal with initial- and final-state continuum wave- 
functions, coefficients no longer have any dependence on principal quantum numbers n. Nor 
do they take on discrete values since they are now functions of the electron's initial energy 
-E, which is continuous. Since we approximate the reaction rate as being strongly peaked at 
the isomer's excitation energy {i.e. Eq. 0}, the reaction rate is simply 

(^i.e.s.^J^m) 1 ^lo^ 

Lastly, we assume the outgoing electron to be in an s-wave state. Table IH] tabulates the 
inelastic scattering coefficients at each isomer's excitation energy Em. Table ITTTl shows rates 
for this reaction channel using the same conditions for the electron fiux as in the previous 
section. 



D. Coulomb excitation via heavy-ion scattering 



Isomeric excitation could also happen via heavy-ion Coulomb excitation. In this case, 
the kinetic energy of the ions is small compared to its rest mass. Hence the ions can be 
treated non-relativistically. Reaction rates for heavy-ion Coulomb excitations can be directly 
calculated using similar methods discussed in the previous sections. However, it is simpler 
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(and perhaps more enlightening) to consider the 'adiabaticity' parameter ja |l 

i = ri^. (19) 
where E is the incoming ion energy and rj is the classical Sommerfeld number given by 

i] = aZiZ2-. (20) 

V 

Here Zi and Z2 are the charges of the two scattering ions. The parameter ^ is a measure 
of the ratio between the collision time and the nuclear period. The larger the value of ^, 
the more adiabatic the process becomes. Under conditions where the reaction is highly 
adiabatic, excitation of the nucleus is unlikely due to the fact that the nucleus has time 
to equilibrate under the influence of the external Coulomb field. Assuming the ions are 
in thermal equilibrium, the energy E can be approximated by \kT . Considering like-ion 
scattering, i.e. Z\ = Z2, the range of temperatures considered in this paper give 

10^ < ^ < 10^. 

The cross section for heavy- ion scattering scales roughly as a ~ e~^'^^j3|; which in turn gives 
a similar scaling for the reaction rate. This rate is extremely small for large ^. It is interesting 
to note that the exponential suppression is primarily a property of the repulsive Coulomb 
field jl^. The Coulomb barrier of these ions is too strong at these plasma temperatures 
for any isomeric excitation. For an attractive field with Zi = 1 (as in inelastic electron 
scattering), this exponential behavior is not as pronounced. 



E. N(7,7')N* and N(e,e'7)N* 

Another possible excitation process is through excitation to some high lying state, fol- 
lowed by gamma decay to the isomer. Excitation can be done via all the processes mentioned 
above, and in particular through (7,7') and (e,e'7) reactions. If one assumes that the rate 
of decay to the isomer is fast compared to the excitation process, then the rate of producing 
the isomer is approximately equal to the rate of excitation to the high-lying state. At first 
glance, such a process might seem plausible, since many of these higher-lying states are 
connected to the ground state by transitions of low multipolarity (A <3,4). Such proba- 
bility transitions may be larger by an order of magnitude or so compared to the isomeric 
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transitions. However, the fluence of incoming particles with the necessary energy to allow 
excitation is exponentially suppressed. Hence the gain from going to a lower multipolarity 
transition is quickly lost due to the fact that these high-lying states have large excitation 
energies compared to the isomers. 

For example, consider the case of fgY. The isomer state is connected via Ml transition to 
the 7/2+ state at energy 2529.87 keV. The transition probability between these two states 
is large (ti/2=-08 ps). The 7/2+ state can be connected to the ground state via an M3 
transition. But since the energy of this high lying state is large compared to the isomer, 
the rate of populating this state is much smaller than for the rate of directly populating the 
isomer. This suppression is roughly 

-^7/2+ _ (,~^530/kT ^ ,Q21/kT 
R. ^ p-909/kT 

Here kT is in units of keV and Risomer is the rate for producing the isomer directly from the 
ground state (via M4 transition). For a plasma temperature of 1 keV/10 keV, this suppres- 
sion factor is vanishingly small. Similar suppression factors are found for the other isomers 
considered in this paper. Hence we conclude that population of isomers via (7, 7') or (e, e'7) 
excitation is extremely unlikely (as opposed to the case of neutron-induced excitations, for 
example). 



III. CONCLUSION 



Nuclear isomers have played an increasingly important role in both astrophysics and ap- 
plied physics phenemona. Typically these isomers are assumed to be populated solely by 
strong nuclear reactions. However, if the temperatures of the environment in which the 
isomers exist are high enough, there is the possibility that they be populated electromagnet- 
ically as well. In this paper we addressed this latter issue for typical plasma temperatures 
found in mundane laboratories and astrophysical environments {i.e. kT~ 1 to 10 KeV). 

In Sect. In] we gave expressions for the reaction rates of the following electromagnetic 
processes: photo-absorption, inverse internal conversion, inelastic electron scattering, and 
Coulomb excitation. In Sect. Ill El we argued that rates for (7,7') and (e, e'7) excitation 
were negligible. We made the assumption that the nucleus was a point particle and was 
completely stripped of its electrons due to the plasma environment. Furthermore, we only 
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considered s-wave scattering since this is usually the dominant channel. We found that all 
these processes contributed negligibly to the overall reaction rates for populating the nuclear 
isomers listed in Table UJ Table UTTl summarizes our claims. 

These small rates can also be understood from simple dimensional analysis. As an exam- 
ple, consider photo-absorption on the ground state of ^•^^U. Using Eqns. [71121 as guides, naive 
dimensional analysis gives the reaction rate for producing Uranium isomers as 

^ |: . ^, (21) 

where in this case /=3 and Q = ~ 2.4 x 10~®. The dimensionless parameter Q sets the 
scale of the transition. Similar expressions for internal conversion and electron scattering 
can be determined from dimensional analysis as well, though in these cases it is the electron's 
momentum that determines the size of Q and not the photon's momentum. In all of these 
cases, Q is very small. An exception would be Coulomb excitation. However, as already 



mentioned in Sect. IIIDt the Coulomb barrier is prohibitive at the temperatures considered 
here. 

Clearly the assumption that the nucleus is completely stripped of electrons is not valid. 
Furthermore, the expressions given in Sect. |nl were made under the assumption that the 
fiuence of incoming particles is in thermal equilibrium. This is most likely not true in 
applied physics applications. However, due to the extremely small rates listed in Table UTTl 
more sophisticated calculations that take these issues into consideration are not warranted. 
Hence the assumption that nuclear isomers are produced entirely via strong nuclear reactions 
for the range of temperatures considered in this paper is valid. 
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TABLE I: Nuclear isomers and their properties. 



Nucleus 


ground state (J'^) 


multipole 


isomer (J'^) 


excitation energy 


n/2 


87y 


1/2- 


M4 


9/2+ 


380.7 keV 


13.37 hr 


SSy 


4- 


E3 


1+ 


392.9 keV 


0.3 ms 


SSy 


4- 


M4 


8+ 


674.6 keV 


13.9 ms 


89y 


1/2- 


M4 


9/2+ 


909 keV 


16.06 s 




0+ 


M8 


8- 


1.174 MeV 


4 s 




0+ 


E16 


16+ 


2.446 MeV 


31 yr 


193 


3/2+ 


M4 


11/2" 


80.2 keV 


10.53 d 


235 -y 


7/2- 


E3 


1/2+ 


77 eV 


26 min 


Am 


1- 


E4 


5- 


48.63 keV 


141 yr 
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TABLE II: Internal conversion and inelastic scattering coefficients. The second column gives the 
orbital in which the electrons are captured during inverse internal reactions. The third column lists 
the conversion coefficients calculated from relativistic hydrogenic wavefunctions, (i.e. Eos. llUltT^ . 
The fourth column gives conversion coefficients derived from Dirac-Hartree-Fock calculations. The 
last column gives scattering coefficients for s-wave electrons calculated at the isomer's excitation 
energy Em- 



Isomer 


final electron state 


ai,c. (r.h.wfs) 


ai.c. (D.H.F.)'' 


ai.e.s.(-E'm) 


87mY 


lSl/2 


.227 




3.54 


SSmly 


lSl/2 


.0245 




1.48 


88m2Y 


lSl/2 


.0226 




.323 


89mY 


lSl/2 


.0078 


.0819 


.107 


178ml jj£ 


lSl/2 


.342 




3.05 


178m2jj£ 


lSl/2 


.109 




20.47 


193mjj, 


2Si/2 


3641 


2610 


1.9E5 


235m^ 


39si/2 


1.5E17 


2.62E17'' 


1.9E22 


242m Ana 


2Si/2 


3724 




4.5E5 



^Ovi shell result 
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TABLE III: Order of magnitude rates for electromagnetic reactions involving photoabsorption 
(P.A.), inverse internal conversion (LLC), inelastic electron scattering (I.E.S.), and Coulomb exci- 
tation (C.E.). Rates are in units of inverse seconds, and come in pairs. The first number corresponds 
to the rate at kT = 1 keV, while the second to kT = 10 keV. We do not list rates for (7, 7') and 
(e, e'7) since these are vanishingly small as well (see discussion in Sect. UTE)) . 



Isomer 


P.A. 


LLC. 


I.E.S. 


C.E. 


87mY 


10- 


-171 


10- 


-22 


10-173 


10- 


-26 


10-172 


10 


-25 




e- 


10^1 


e- 


10^ 


88ml Y 


10- 


-168 


10- 


-14 


10-172 


10- 


-20 


10-170 


10 


-18 




e- 


10^1 


e- 


10^ 


88m2Y 


10- 


-296 


10- 


-32 


10-300 


10- 


-38 


10-299 


10- 


-37 




e- 


10^1 


e- 


10^ 


89mY 


10- 


-397 


10- 


-41 


10-401 


10" 


-47 


10-399 


10- 


-46 




e~ 


10^1 


e~ 


10'^ 


178ml jj£ 


10- 


-505 


10- 


-70 


10-522 


10- 


-75 


10-521 


10- 


-74 




e~ 


108 1 


e~ 


10^ 


178m2jj£ 


10" 


1116 


10- 


-160 


10-1119 


10- 


-165 


10-1118 


10" 


-164 




e~ 


10^1 


e~ 


10"^ 


193m jj. 


10 


-46 


10" 


14 


10-44 1 


10- 


14 


10-42 1 


10- 


12 




e- 


10^1 


e- 


10^ 


235m^ 


10 


-21 


10- 


20 


10-9| 


10- 


7 


io-4| 


10- 


2 




e- 


10*1 


e- 


10^ 


242"^ Am 


10 


-31 


10- 


12 


10-30 1 


10- 


12 


10-28 1 


10- 


10 




e' 


10^1 


e~ 


10'^ 
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(e) (f) 

FIG. 1: Feynman diagrams representing the electromagnetic processes considered in this paper. 
Plot (a) represents photo-excitation, (b) represents inverse internal conversion, where the final 
electron is bound, (c) shows electron scattering, where the final electron belongs to the continuum, 
(d) shows ion-excitation, where N' is another ion. (e) shows (7,7') excitation. Finally, (f) shows 
(e, e'7) excitation. 
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